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τε, INTRODUCTION 


A. WHY MESOSCALE MODELS ARE IMPORTANT 

Scientists and weather forecasters are often 
interested in understanding small-scale phenomena within 
specific regional areas. In many cases, intended areas of 
study would be impossible without the employment of 
mesoscale models.  Mesoscale models take accepted global- 
scale atmospheric models and, with more detailed physics 
and smaller grid spacing, provide details of smaller scale 
atmospheric phenomena smoothed over by the global-scale 
models. This ability to better resolve topography and look 
with finer spatial and temporal resolution makes a 
mesoscale model an invaluable tool in many places. For 
example, in areas of complex terrain, such as the Los 
Angeles Basin, accurate representation of wind flow 
requires a scale of 4-5 km, or less (Stauffer et al., 
2000). This requirement far exceeds the resolution of any 
Global scale. Without mesoscale models, better 
understanding of atmospheric structures within some areas 
of study would be impossible. 

For the U. S. Navy, mesoscale models provide valuable 
and necessary information into Tactical Decision Aids 
(TDAs). Whereas atmospheric soundings can be input to a 
TDA to gain a real time analysis of the atmospheric 
environment, output from mesoscale models is used to gain 
insight of the atmosphere in the future. Mesoscale models 
are more likely to capture atmospheric details that a 
rawindsonde would show than a synoptic-scale model would. 
Assuming the mesoscale model has been initialized with 


accurate synoptic-scale structure, the mesoscale model 


1 


output gives the warfighter the opportunity to better 
exploit the environmental battlespace in the short-term 


future. 


В WHY FORECASTING AIR QUALITY IS IMPORTANT 

The California Air Resources Board (CARB) was a 
cosponsoring agency of the South Coast Ozone Study of 1997 
(SCOS97). One objective of CARB and its cosponsors in 
SCOS97 was to measure particulate levels over the Los 
Angeles Basin. Many local, regional, and national agencies 
are interested in the measurement of atmospheric 
particulates and the short-term prediction of the 
atmospheric pollutant contents. Locally, hospitals may be 
concerned with poor air quality days because it could mean 
increased numbers of patients with respiratory problems 
exacerbated by high ozone levels. Regional authorities 
would be interested in high ozone level predictions. With 
high ozone predictions in the short term, appropriate 
officials would issue carpool advisories or public notices 
of increased ozone level risk. By using a pollutant level 
prediction, high concentrations of ozone may be mitigated 
before they actually happen. National and international 
agencies, such as the North American Research Strategy for 
Tropospheric Ozone (NARSTO) have vested interest in ozone 
measurements and predictions.  Forecasted areas of high 
ozone concentration would lead them to better place 
measurement instruments. The U.S. Navy also cosponsored 
the SCOS97 study. Tracking pollutants in the Los Angeles 
basin has direct correlations to tracking smoke, 


biological, or chemical parcels imma liteeral™ satevespaces 


Air quality models (AQMs) are necessary for these 
applications and externally driven meteorological models 
influence these AQMs. The Pennsylvania State University 
(PSU) /National Center for Atmospheric Research (NCAR) 
Mesoscale Model Version 5 (MM5) is recognized as a superior 
model for driving other AQMs and air quality applications 


(Pielke et al., 1991). 


er PREDICTABILITY THEORY 

Atmospheric predictability is the ability of a 
numerical model to properly describe the atmosphere at a 
future time. All models, in an attempt to maximize 
predictability, strive for accuracy and skill. Accuracy is 
an objective measure of how well a particular model’s 
forecast fits the observed sequence of events. Skill is a 
relative measure of forecast accuracy where one particular 
forecast is compared to another forecast (Nuss, 2000). 

Figure 1-1 illustrates the hypothetical error of a 
numerical model forecast compared against a climatology- 
based or persistence-based forecast. A climatology-based 
forecast has constant error over time. A persistence-based 
forecast has less error initially until its error degrades 
to that of climatology at time t and is worse than 
climatology beyond t2. The numerical model forecast has 
relatively larger error in the very short term. This 
initial error is due to the lack of consistency between the 
initial model state and the dynamic allowances set in the 
model physics (Kuypers, 2000). Initial error is also due 
to unresolvable scales of motion the model cannot capture 
(Anthes and Baumhefner, 1984). After the model adjusts and 


achieves a stable agreement between the initial 
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observations and model physics (also called "spin up"), the 
process leading to the error peak (at ti) decreases to a 
minimum. Here the model shows skill over both climatology 
and persistence. The initial errors grow over time, and at 
some point t3, the model no longer shows skill over 
climatology. 

Many factors affect atmospheric predictability, one of 
the most important being resolution. Mesoscale models with 
finer resolution grasp more spatial and topographic detail, 
and, using more sophisticated physics, can incorporate 
smaller scales of motion that a larger scale model could 
not resolve. When applied over a tightly defined 
geographic region, the mesoscale model can resolve motions 
caused by local forcing, necessary for mesoscale phenomena 
that a global or coarser model would not resolve. Over 
complex terrain, for example, a model needs a resolution of 
4-5 km or less in order to properly resolve the flow 
(Stauffer et al., 2000). 

For a given area of study, if weather is well 
observed, we expect less initial error with a finer 
resolution model. Finer resolution is more likely to 
capture smaller-scale phenomena, and a finer resolution 
model will maintain these phenomena longer through a 
forecast period. However, with more detailed physics and 
finer spatial separation, mesoscale models can have great 
variability over an area. 

These phenomena of interest have shorter lifespans 
than those of synoptic-scale models, adding more 
complications to mesoscale forecasts. It is difficult to 
observe these phenomena on such fine spatial- and time- 


scales, let alone forecast them. Observations aloft are 
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even sparser than surface observation sites, making the 
ability to capture accurate vertical structure very 
challenging. 

It is fair to compare one model's forecasts under 
different resolutions provided that a reasonable time scale 
is considered for all of the forecasts. Error growth, at 
any resolution, eventually contaminates all scales of 
motion (Anthes and Baumhefner, 1984). To minimize this, we 
consider only short-term forecasts (specifically 3hr, 6hr, 


Shu and 112г) femthis study. 


D. MODEL VERIFICATION 

Model verification is most often done one of two ways: 
comparing model forecasts against model-generated analyses 
or comparing model forecasts against observations. The 
former is most convenient because the forecast and analysis 
fields are on collocated model grids. The latter is more 
challenging because there are rarely enough observations 
over a particular area to create a picture as complete as 
the model-generated analysis. In areas of sparse 
observational data, the first guess forecast is the 
heaviest influence on the current analysis. The Cooperative 
Program for Meteorology, Education and Training (COMET) 
Numerical Weather Prediction Module discusses the details 
of each model verification method. 

This case study compares synoptic patterns from NOGAPS 
(Navy Operational Global Atmospheric Prediction System) to 
MM5 to ensure MM5 has the correct "big picture" (not 
included) and compares model forecasts against observations 
to generate statistical results. The SCOS97 provides an 


exceptionally rich data field over the area of study. 
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Е. HYPOTHESIS 

Meteorologists generally accept that a static 
initialization, or cold start, of a mesoscale model 
produces poorer forecasts than a model with a dynamic 
initialization, or warm start. A cold start of a mesoscale 
model incorporates a global model (in this study NOGAPS 
1%x1 horizontal resolution) and National Weather Service 
(NWS) observations for defining initial and lateral 
boundary conditions (LBCs). From these inputs, the 
mesoscale model derives its analysis and forecasts. A warm 
start uses the same, but also incorporates a previous 
mesoscale model forecast as its first guess. Conventional 
wisdom says the warm start is better because it maintains 
the mesoscale structure generated by the first guess that a 
global model cannot. Essentially, the cold start resets 
the mesoscale structure and must redevelop it at each 
initialization. 

Also, a model with more complex physics and complex 
terrain resolution should perform better than a simpler 
model. Sophisticated physics and terrain resolution should 
allow a model to see and interpret finer mesoscale 
phenomena. It follows that the more sophisticated model 
physics, for an accurately defined initial state, should 
have less initial error. Therefore, a model run with 
sophisticated physics has less initial error propagating 
through its forecast cycle. Also, the model may trace a 
parcel more precisely over terrain through a forecast, 
lessening the model user’s need to round off or interpret 
some values. 

Another widely accepted theory is more observations 


fed into a model at analysis time will lead to better 
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forecasts. Again, more observations lead to a more 
complete picture of the atmosphere at analysis time. Any 
model analysis will have some initial error. However, 
decreasing initial error means less initial error will 
propagate, yielding better forecasts. 

We hypothesize these principles are not always true. 
Because the phenomena a mesoscale model attempts to follow 
are short lived in both time and space, there will always 
be error in any analysis or forecast. А warm start may 
propagate an error in a model's forecast further than would 
a cold start. More sophisticated physics may not provide 
any significant benefit to a forecast if the initial error 
is too great or if the model physics use assumptions not 
valid for the scale of interest. Lastly, more observations 
may not add any significant insight to the initial state of 
the atmosphere the model is trying to resolve because they 
might not be located in a position relevant to defining 


important small-scale structural details. 


F. OBJECTIVES 

In this study, we examine four one-way nested domains 
of MM5, namely MM5 81km, 27km, 9km, and 3km. We examine 
the ability of each domain to capture the flow over the 
area of study and examine the error statistics of each 
model compared against observations. We will run a control 
experiment model run for each domain using a warm start and 
sophisticated planetary boundary layer (PBL) scheme. Our 
first experiment will be to run a cold start for each 
domain, maintaining the same complex PBL scheme. Our 
second experiment will look at running a warm start using 


simpler physics (a less complex PBL scheme). Lastly, we 
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will run the same control experiment and inject extra 
observations from the SCOS97 study into the model 
initialization that were not part of the original control 
case. 

We will show that, for our case, the accepted 
conventional wisdom mentioned earlier is not necessarily 
correctas 

A description of the model used in this study will be 
given in Section II. Background information related to the 
SCOS97 study and related observations will be given in 
Section III. Section IV will discuss weather over the area 
of study for 4-7 August 1997. Our methodology and results 
will be discussed in Section V and Section VI. Lastly, 
Section VII will offer our conclusions and ideas for 


further study. 


ТТ. MODEL DESCRIPTION 


The PSU/ NCAR MM5 Version 3.3 is used for this study. 
MM5 is a one or two-way continuous, multi-nested, non- 
hydrostatic, terrain-following, sigma coordinate, primitive 
equation model. For this study, MM5 was run on a one-way 
continuous, quadruply-nested configuration with horizontal 
grid spacings of 81km, 27km, 9km and 3km. The 81km domain 
is centered at 34.2N 119.38W with 41 points in the E-W 
direction and 45 points in the N-S direction. The 27km 
domain has dimensions 43 x 43. The 9km domain has 
dimensions 73 x 73 and the 3km has dimensions 97 x 97. 
Figure 2-1 shows the horizontal extent of each domain. 

NOGAPS with 1? x 1? resolution is used for both the 
initial and boundary conditions. Warner and his coauthors 
describe limitations and concerns about lateral boundary 
conditions applied to local area models (Warner et al., 
1997). 

Two-dimensional multiquadric interpolation (2DMQ) 
(Nuss and Titley, 1994) is used to blend first guess grid 
fields and observations to MM5 grid fields. There are 
several key points to mention about the 2DMQ method. 
First, as the name implies, the method is two-dimensional. 
Each isobaric level is solved independently of any other. 
For example, wind direction at 850mb is solved for 
independently of wind direction at 925mb. 

A second point to note is this method is univariate, 
meaning each parameter in the model is solved for 
independently of another. For another example, wind speed 


is solved independently of temperature. 


One compelling argument for using the MM5 is the 
ability to test the sensitivity of model forecasts to 
different physical parameterization schemes. We choose to 
compare the more sophisticated physics in the Gayno-Seaman 
PBL scheme (Gayno et al., 1994) to the simpler physics of 
the Blackadar PBL scheme (Blackadar, 1979). We use a five- 
layer soil model parameterization scheme and no convective 
parameterization scheme for our MM5 runs. 

All domains incorporated a vertical resolution of 30 
vertical levels. A total of 12 levels are located between 
the 850mb level and the surface. 

Topographic resolution for each domain is incorporated 
from the U.S. Geological Survey. The 81km and 27km domains 
derive their topography and land use from 19km terrain and 
land use resolution. The 9km domain derives its topography 
and land use from 9km terrain and land use resolution, and 
the 3km domain derives its topography and land use from 1km 


terrain and land use resolution. 
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ТТТ. BACKGROUND 


А. ABOUT THE SCOS97 STUDY 

SCOS97 incorporated a very large base of 
meteorological observations spread over an 80,000 square 
mile area within the Los Angeles basin. The massive amount 
of data was collected over a period from 16JUN1997 to 
150CT1997 from 259 different observation sites. 

This paper focuses on the period from 04-07 August. 
These dates envelop a measured high ozone event, with 
surface ozone measurements peaking near 2200Z on August gth 


(Boucouvala, 2000). 


EL DESCRIPTION OF OBSERVATION STATIONS 

Figure 3-1 depicts the general topography of the Los 
Angeles basin. The basin is bordered on the west and south 
by the Pacific Ocean. The Santa Monica and San Gabriel 
Mountains border the northern area of the LA basin and the 
San Bernadino and San Jacinto mountains mark the eastern 
boundary. The Santa Ana Mountains are southeast of Los 
Angeles and also affect area weather. Figure 3-2 shows the 
boundaries of our area of verification and also shows how 
each of the 37 chosen observation stations fits with the 
local terrain. Table 3-1 provides the location 
specifications of each station. 

Surface sites provided hourly observations. Five 
rawindsonde sites made observations every six hours (002, 
062, 122, 182). Continuous profiler data (mean values 
only) came from nine sites. All observations were quality 
checked by CARB, and all observations of uncertain quality 
were excluded from our study. 


ШЕ 


All observations were noted at heights in meters above 
ground level (AGL). These values were used as the 
“control” for the model measurements. Since MM5 is a sigma 
coordinate model, all model values were interpolated 
vertically to the same level as the observations for 
comparison purposes, with the exception of the surface 
observation comparison. Surface comparison was done using 
the model forecast parameters at the lowest sigma level, 
located at about 22m AGL. 

Infrequent gaps in the observation data exist for 
varıous reasons. The most likely reason for these gaps is 
due to the quality control scrub by CARB. We should also 
note the gaps in observational data at non-synoptic times. 
At the synoptic times, we are afforded the additional data 


from the rawindsonde measurements. 
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ТУ. METHODS 


А. HOW WE SEGMENTED OUR DATA 

Table 3-1 and Figure 3-2 describe the stations from 
which our observational data came. We divided our 
atmospheric area of interest into six layers, specifically 
described in Table 4-1. Over the area of study and within 
each layer, we summed all appropriate data points to derive 
an average for a particular time. For any particular time, 
there will be six values for a given parameter, one value 
for each of the six layers. 

The same principles were applied to model forecast 
results, and the MM5 data was interpolated to the exact AGL 
location as the upper air observation data. The 
interpolated values come from model grid point values 
within the verification area. The 81km domain has three 
grid points over this area, while the 27km, 9km, and 3km 
domains have 36, 322, and 2912 grid points, respectively. 

This study is concerned with short-term forecasts, 
namely MM5's 3hr, 6hr, 9hr, and 12hr forecasts. These are 
the only forecasts plotted over time. To further clarify, 
examine Figure 4-1. At 05/032 апа 05/152, five values for 
mean temperature in the surface layer are plotted, namely 
the observed mean temperature and the mean temperature 3hr 
forecasts for all four domains. At 05/06 and 05/18Z, the 
mean observed temperature is plotted along with the 6hr 
forecasts of the four MM5 domains. Likewise, 05/09Z and 
05/212 plot 9hr forecasts, and 05/12Z апа 06/002 plot 12hr 
forecasts. 

From the observational data, we compute mean and 


standard deviation fields for wind speed, wind direction, 
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temperature, and virtual temperature and derive the same 
information from the model forecasts. In addition, we 
calculated a root mean square difference and bias for each 
model domain. We used the same definitions of error and 
bias as the U. S. Navy's Fleet Numerical Meteorology and 
Oceanography Center (FNMOC). Error is the difference of 
model forecast minus observation. Bias is model's tendency 
to over or under-forecast a parameter.  FNMOC's Monthly 
Models Report has a section of routine statistical measures 
that describes these in greater detail. 

A parallel study (Mohammed, 2000) uses the same data 
set to study the WOCSS (Winds Over Critical Streamline 
Surfaces) model. Some data related to his study is also 


seen in this study's results. 


P. WARM START VS. COLD START 

Our control experiment (noted CTRL), the normal ММ5 
set up to which we compare our other experiments, is a warm 
start. It is initialized using a previous 12hr ММ5 
forecast blended with NWS surface and rawinsonde 
observations. The MM5 12hr pre-forecast is generated using 
a cold start initialized with NOGAPS analyses and available 
NWS observations. This CTRL run employs the sophisticated 
Gayno-Seaman PBL scheme (Gayno et al., 1994). 

The first study is a comparison to a cold start of MM5 
(noted COLD). In COLD, MM5 is initialized using only 
NOGAPS analyses and NWS observations. COLD also employs 


the Gayno-Seaman PBL scheme. 
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Gu VARYING PBL SCHEMES 

The second experiment (noted TEST1) looks at changing 
the PBL scheme used in MM5. For this study, we ran MM5 
using the Blackadar PBL scheme (Blackadar, 1979). The 
Blackadar PBL scheme is a lower order closure model than 
the Gayno-Seaman PBL scheme. All other aspects of the 
model set-up in TEST1 match those in CTRL. 


р. INJECTING ADDITIONAL OBSERVATION DATA 

The final experiment (noted TEST2) examines the effect 
of injecting additional observations into MM5’s 
initialization. The NWS observations are the normal input 
expected to initialize MM5 on a regular, consistent basis. 
Our study takes additional observations from the SCOS97 
data set and inputs them into MM5's initialization. These 
observations were not part of the input to CTRL and are 


independent from those used in model verification. 


E. DATA PROCESSING 

Weather chart analysis was done using GARP 
Version 5.4.  FORTRAN programs were used to compute some 
observation and model statistics and to bin the data into 
appropriate files. Microsoft Excel computed additional 


statistics and was used for plotting statistical results. 


US 


THIS PAGE INTENTIONALLY LEFT BLANK 


15 


У. WEATHER DISCUSSION 


A. CLIMATOLOGICAL DISCUSSION OF THE LOS ANGELES BASIN 

The Los Angeles basin has a climate dominated by two 
distinct seasonal rainfall patterns, a rainy winter season 
and a dry summer season. The modified Kóppen 
classification system categorizes this climate as 
Mediterranean (Bruno, 2000). This study is embedded in the 
heart of the dry summer season. During this season, 
southern California is dominated synoptically by the 
Eastern Pacific High. Warmer, dryer air aloft subsides 
over cooler, moister air at the surface to create a marine 
layer, the predominant weather feature during the Los 
Angeles summer. Under this high pressure, synoptic winds 
are generally light and variable. Mostly clear skies over 
inland areas and often times cloudy skies over coastal 
areas lead to a substantial surface temperature gradient. 
This temperature front induces the prevailing summertime 
sea/land breeze regime over the Los Angeles basin. 

Figure 5-1 depicts a typical diurnal sea breeze regime 
(Ahrens, 1991). Uneven heating rates between land and 
water induce a mesoscale circulation, producing onshore 
winds during the day (sea breeze) and offshore winds during 
the night (land breeze). An important consequence of the 
sea breeze is the sea breeze front, the onshore push of 
cooler marine air over land. Ahrens provides more details 
of the sea/land breeze regime in his text. 

The rainy winter season (NOV-APR) is influenced by the 
migratoryetrack of the Polar’ Front-Jet.s Cool air ο ας 
storms in the mid-latitude cyclone track break down any 


substantial marine layer. This marine layer is a concern 


ШЕГІ 


because it tends Бомашар ап εποσ. It logically 
Follows that more bad air quality days exist in the summer 
months than in the winter months. The National Weather 
Service publication NWS WR-261 is an excellent data source 
for providing much greater detail on the climatology of Los 


Angeles. 


В» SYNOPTIC SITUATION 

For the period between 040000ZAUG1997 and 
080000ZAUG1997, large scale synoptic conditions are 
generally unchanging. For readability, all ddhhmmZMONYEAR 
times shall be abbreviated. For example, 051200ZAUG1997 
shall be noted as 05/122. 

Building heights at 300mb over the western U.S. keeps 
the long wave trough and the Polar Front Jet far from any 
areas that may affect southern California. In fact, NOGAPS 
analyses show no 300mb winds greater than 15-20kts over 
southern California during this period (not shown). 

АС 04/0027, a 5580m cut-off low sits west of Vancouver 
Island and a strong ridge sits over the western U.S. at 
500mb (Figure 5-2). The low fills slightly over the next 
36 hours (Figures 5-3 through 5-5) while heights rise over 
southern California and the western U.S. At 06/002, an 
approaching low over the Aleutian Chain influences the 
Vancouver low (Figure 5-6) while the western U.S. ridge 
reaches its maximum heights. Тһе long wave trough deepens 
south of the Aleutian chain through 08/0027 while the ridge 
over California flattens and slightly weakens (Figures 5-7 
through 5-10). Weak and constant vorticity values over 
southern California during this period show little, if any, 
vertical motion exists at this level. 
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Тре 925mb series of geopotential height and 
temperature plots shows a general warming trend throughout 
the period of interest.  04/00Z shows a closed 27C isotherm 
over southern California (Figure 5-11), and the chart 
following 12 hours later (Figure 5-12) shows that pool of 
warm air 3C cooler and closer the California/Mexico border. 
The charts on August 5* show not only a broader span of 
warmer temperatures over California, but also closed 
isotherms 3C greater over southern California on both the 
002 апа 12Z analysis. 05/002 shows another 3C change over 
California near Point Conception, as well as a 3C change 
over the California/Arizona border, increasing to closed 
isotherms of 33C and 36C respectively (Figure 5-13). The 
122 analysis on August 5th (Figure 5-14) shows a 3C warm up 
over southern California bight region from the previous 24 
hour 12Z time frame. Temperatures at 00Z on August ens 
(Figure 5-15) remain fairly constant compared with the 
previous day, while 122 temperatures (Figure 5-16) over the 
bight region cool slightly from the previous evening, 
evidenced by the lack of a closed 30C closed isotherm. 
Temperatures over the bight region cool by a few degrees on 
the 07/00Z analysis (Figure 5-17), while the warm 
temperatures over the California/ Arizona border change 
very little from the previous day, hinting toward the end 
of the warming trend. Temperatures at 07/122 (Figure 5-18) 
Similar to those from the previous day further support the 
conclusion of the warming trend. Temperatures continue to 
cool at 08/00Z (Figure 5-19), with the warmest pocket of 


warm air shifting eastward into Arizona, drawing the 


thermal tongue eastward and allowing temperatures over the 
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bight region to cool another 2-3C from the previous 24 hour 
period. 

The surface plot of mean sea level pressure shows the 
diurnal shift of the inverted thermal trough over 
California. Figures 5-20 through 5-28 show the series of 
surface pressure analyses. Intense heating during the day 
draws the thermal trough further east from the coast, 
evidenced by the 00272 plots. At night (12Z plots), the 
trough is slightly less pronounced and further toward the 
Pacific coast. Surface pressure change over southern 
California is minimal during the period. However, the 
eastern Pacific surface high pressure center west of 
central California along 130W longitude builds from 1024mb 
оп 05/007 το 10261 оп 06/007. This time corresponds to 
the peak of the upper level ridge as well as the lowest 
surface pressure over California/ Arizona/ Mexico border. 
Ihe result is an increased pressure gradient over the 
central California coast for a brief time between 00Z and 
12Z on August 6° (Figures 5-24 and 5-25). The surface high 
continues to build through the rest of the period, but 
shifts northward and westward as it does so, producing 
little impact on any pressure gradient over central or 


southern California after 06/122. 


Co MESOSCALE SITUATION 

For the mesoscale meteorological situation, attention 
is focused on the period from 04/122 to 07/002 

We make some assumptions and define thresholds for 
winds over the area of study. We assume strong onshore 
flow (sea breeze) to be a primarily west-southwesterly 


wind, from 220° to 290°, and offshore flow (land breeze) to 
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be a primarily north-northeasterly wind, from 020? to 100?. 


1 


More importantly, we define a speed of 3 ms or less to be 


an insignificant enough wind speed to determine a 
verifiable wind direction. That is, 3 ms! winds are 
considered, for this study, light and variable. 

The plots of domain-averaged mean observed 
temperatures across the timeline show a distinct difference 
in between the lower layers and the upper layers (Figure 5- 
29). The surface layer demonstrates a clear diurnal cycle, 
warmer temperatures at the 18 and 0027 time frames and 
cooler temperatures at 0627 апа 122. Layers two and three 
also show a diurnal temperature wave, and though these two 
layers are generally in phase, they lag the surface plot. 
Further aloft, the recognizable diurnal temperature wave 
becomes a bit more complex. 

All layers seem to verify the general warming trend 
seen on the synoptic charts. Over the area of study, the 
general warming trend begins on August 4*5 and continues 
until early morning of Augustwe:^ wea Ἱ layers show a warmer 
temperature at 05/12Z than the time 24 hours previous. The 
same is true when looking at the warm part of the diurnal 
cycle. Five of the six layers at 06/00Z are warmer than 
the 05/0027 piloti 

After 06/12Z, the warming pattern of the 48 hours 
earlier starts to disappea At о и Нонна неее ыы | 
surface layer and layer four have continued to warm up 
compared to the previous evening, but the other four layers 
have returned to their previous low or have cooled. 

Looking at wind speeds for the time period (Figure 5- 
30), surface wind speeds peak at approximately 21Z of each 


day and slow to a minimum at night between 0927 апа 122. 
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Peak upper layer winds show a more complex variation from 
diurnal change. Speeds on August 4" апа 5° аге relatively 
weak, rarely getting aboves5 ms in any layer at any time. 

For August 4°", the peak time of the sea breeze occurs 
between 21Z and 002. Looking at the 05/00Z time frame, 
winds are strongly onshore from the surface up through 
layer four. Layers five and six show light and variable 
winds. Recalling Figure 5-1, we conclude that our six 
layers fail to capture the entire sea breeze solenoid, as 
no strong evidence of return flow exists (Figure 5-31). 

At 05/12Z, winds are light and variable throughout all 
six layers. Like the night before, no evidence exists of 
any organized land breeze. 

Wind speeds increase again as temperatures rise, and 
the sea breeze again dominates the surface layer by 05/212. 
At 06/00Z, wind speeds in the upper layers follow the 
surface and strengthen, but they are greater than the winds 
from the previous afternoon. In fact, layers three through 
six increase by more than 2ms'' from their wind speed on the 
05/002 chart. Conversely, mean speeds in the surface layer 
and layer two decrease. 

Comparing the 06/00Z observations against those from 
24 hours previous reveals a definitive change in the 
strength of the solenoid. Peak temperatures at 06/002 
would drive a stronger sea breeze in the surface layer, not 
a weaker one. Upper layers would also tend to reflect a 
mean direction similar to the sea breeze of the day before. 
Instead, upper layer winds have gained a more northerly 
component at 06/00Z (Figure 5-31). A decrease in wind 
direction standard deviation over layers three through six 


(not shown) indicates the wind is more consistently from 
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one direction over the area of study. It is the brief 
synoptic change interfering with Los Angeles mesoscale 
regime that alters the state of the solenoid early on 
August 6°. This interruption is very brief. Surface winds 
are light and variable again by 06/12Z, and upper level 
winds are decreasing as well. As for the previous night, 
no strong demonstration of a nighttime land breeze exists 
in any layer. 

By 07/00Z, temperatures at the surface begin to 
increase slightly. The surface layer sea breeze strength 
increases, and winds throughout all six layers demonstrate 
onshore flow. 

Another important feature to trace in the mesoscale 
regime is the sea breeze front. The push of the sea breeze 
advects cooler marine air inland, lifting the inversion 
base. We trace the sea breeze front from hourly profiler 
observations, specifically examining profiles taken at 00Z 
and 122. 

Profiler stations not located in close proximity to 
the coast (all but PHE and LAX) show surface-based 
inversions at night (122) and various changes in the 
inversion base height during the day (002). A 002 station 
profiler showing an elevated inversion base as well as 
surface temperatures Similar to those at 122 tells us the 
sea breeze front has reached that station, perhaps even 
pushing further inland. Refer to Figure 3-2 to see the 
location of each profiler site. 

At 04/12Z, only PHE and LAS have inversion bases above 
the surface, showing they are the only stations embedded in 
marine air. АЕ 05/002, the sea breeze front has reached 


USC, SMI, LAX, and TIN, evidenced by a rise in inversion 
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base height and/or a drop in surface temperature. The 
front retreats toward the coast again at night, with only 
LAX showing any signs of marine air at 05/12Z. АЕ 06/002, 
the front has again moved inland as far as EMT, SMI, USC, 
апа ТТМ.  06/12Z2 profilers show the front has again 
retreated back to the coast, with only PHE showing evidence 
of marine air. A key station in these profilers is VNS. 
In this entire time frame, VNS shows neither an elevated 
inversion nor a drop in surface temperature. We conclude 
the sea breeze front never penetrates as far inland as VNS. 
In summary, the sea breeze solenoid is the dominant 
mesoscale feature over Los Angeles for the period, with the 
exception of the brief synoptic influence shortly after 


016.002. 
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УТ. RESULTS 


A. PARAMETER AVERAGE TIME SERIES 

The surface layer presents the greatest variability of 
the six layers, not only between each domain, but also 
between each experiment. Figures 6-1 through 6-3 show the 
time series plots for domain-averaged mean wind speed, mean 
wind direction, and mean temperature in the surface layer 
for CTRL. The same is done for each parameter in COLD 
(Figures 6-4, 6-5, 6-6), TEST1 (Figures 6-7, 6-8, 6-9), and 
TEST2 (Figures 6-10,6-11, 6-12). Mean observations are 
also plotted in each figure and are a good comparison 
marker. 

Looking at the wind speed time series, observations 
show a clear diurnal cycle with one peak speed at the hot 
part of the day. No domain in any of the four experiments 
grasps a similar pattern. Each domain in each experiment 
shows over-forecasted speeds in the cold part of the day. 

Wind direction also has a distinct diurnal wave in the 
observations. Very low wind speeds (less than 1.5 ms!) at 
night make a definitive wind direction difficult to 
ascertain. During the daytime, however, it is clear that a 
sea breeze exists. Each model experiment shows a diurnal 
pattern, although there are some significant differences. 
Each experiment shows the different domains forecasting 
nighttime winds with an offshore component greater than the 
observations. In each experiment, the 3km and 9km 
forecasted directions are more similar to the observations 
and have a more damped wind direction wave than 27km or 


8ikm. 
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The temperature time series present some interesting 
results. As the observations show, there is a strong 
diurnal temperature wave that experiences a gradual warming 
over the period, supporting the trend seen in the NOGAPS 
charts. In each experiment, all domains grasp this diurnal 
pattern, although each domain in each experiment lags the 
observations to some extent. By lag, we mean the tendency 
for the models to forecast a minimum temperature 
approximately three hours later than the observed minimum 
is seen. All domains do an adequate job in forecasting the 
warm part of the day, but all experiments fail to cool the 
surface temperature enough at night. At these cool times 
(062, 092, 122), COLD shows forecasted temperatures 1-4C 
closer to the observed temperature than CTRL, TEST1, and 


TEST2 in each domain. 


B. STUDY PERIOD STATISTICS 

Tables 6-1a, 6-1b, and 6-1c show a further simplified 
representation of CTRL's forecasts in Figures 6-1, 6-2, and 
6-3. The mean and standard deviation for wind speed, wind 
direction, and temperature for each domain's forecast have 
been averaged over the entire timeline. The same was done 
for the observed values. In addition, we calculated error 
(RMSD) and bias for each forecast. With over 500 forecast 
points to calculate each error or bias, we know even minor 
differences between domains are significant. Similar 
calculations have been done for COLD (Tables 6-2a, 6-2b, 6- 
2c), TEST1 (Tables 6-3a, 6-3b, 6-3c), and TEST2 (6-4a, 
6-4b, 6-4c). 

Recall the 81km domain had only three grid points over 


the verification area, meaning that all 500 points for the 
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8lkm domain were interpolated from тапу fewer grid points 
than the other domains. It is unlikely that such coarse 
resolution would capture any significant mesoscale 
structure we are interested in. If this domain appears to 
perform better than other domains, we conclude it is 
getting better error statistics for the wrong reasons. 

Beginning with CTRL, forecasted wind speed gains 
significant improvement at 9km over 27km, and slight 
improvement at 3km over 9km.  Forecasted mean, error 
(RMSD), and bias all improve with finer resolution, 
however, all domains have error greater than the standard 
deviation of observed wind speeds, indicating the forecasts 
show poor forecast skill (Kuypers, 2000). 

The 3km wind direction is the closest to the observed 
mean direction. Each domain shows similar error, all of 
which are greater than the standard deviation of the 
observations. 

The 27km, 9km, and 3km domains have a mean temperature 
more than 3C greater than the observed mean. Again, all 
domains have error greater than the observation's standard 
deviation. 

In COLD's statistics, mean wind speeds for 27km and 
9km domains are closer to the observed mean than in CTRL. 
Error and bias have also improved for these domains. The 
3km wind speed statistics are nearly identical for COLD and 
CLR Ga 

Mean wind direction for COLD in each domain has 
changed more than 20° closer to the observed mean direction. 
Direction error in each domain is also less than CTRL, yet 
none are less than the observed direction standard 
deviation. 
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Mean temperature forecasts have also dropped nearly 3C 
for each domain. This measure, as well as error and bias, 
are much better than CTRL. Each domain's temperature 
forecast error exceeds observational standard deviation. 

In TEST1, the simpler PBL scheme, results for each 
domain's wind speed and direction are slightly poorer than 
CTRL. Within TEST1, 3km and 9km domains show better speed 
and direction statistics than 27km. 

TEST1 temperature statistics show slight improvement 
over CTRL, most noticeably in each domain’s bias. Again, 
3km and 9km prove better than 27km. 

Lastly, TEST2’s forecasts present numbers very similar 
to the statistics of CTRL. Adding extra observational data 
into MM5’s initialization did little to change any values 


seen in the statistical result tables. 


Ge VERTICAL PROFILES 

Figure 6-13 depicts the 9km domain-averaged vertical 
structure for each experiment at 04/12Z. Mean values for 
each parameter were derived at each layer (as done in the 
previous section) and plotted at the mid-level of each 
layer. This was done for each synoptic time through 07/002 
(Figures 6-14 though 6-18), and the same was done for the 
3km domain (Figures 6-19 through 6-23). 

As seen in the earlier times series plots, all domain 
forecasts in each experiment had difficulty in the cool 
part of the day (122). Looking at these times over the 
period of study shows very interesting results, such as the 
strong similarity between the 3km and 9km domains in each 


experiment. 


28 


Focusing near the surface on the 9km 04/12Z profile 
(Figure 6-13), each experiment forecasts temperatures much 
warmer than actually seen in the observations. The next 
two nichts (Figures 6-15 and 6-17) show mean observation 
temperatures increasing, once again illustrating the 
warming trend. Note the model forecasts do not 
significantly change. The model forecasts are getting 
closer to the truth, but it is because the truth is moving 
toward the model solution. Each night, in both the 9km 
domain and 3km domain, the COLD profile comes closest to 
the observed vertical temperature structure. 

Nighttime observations show a temperature inversion 
existing up into the third layer, while both domain 
forecasts only show the inversion existing up into layer 
two. COLD’s profiles in both domains at 05/12Z and 06/122 
are the only profiles that show a PBL height within 500m of 
the observed PBL height. 

Vertical wind speed structure at night shows the same 
forecast problem. For the 9km domain, each experiment 
shows wind speeds which are over-forecast at night. The 
04/122 profile of COLD's wind speed is the only exception. 

Another feature of the profiles is how CTRL, TESTI, 
and TEST2 all show coast-parallel winds in layers three 
through six. Observations show winds with an onshore 
component at these heights AGL. COLD’s wind direction 
profile stands out from the others. 

Considering the water temperature of the Pacific Ocean 
to be essentially constant at 04/12Z, 05/122, апа 06/122, 
the forecast from each domain and experiment shows a much 


stronger coastal temperature gradient than the observed 
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temperatures suggest. This implies a stronger thermal wind 
in the forecasts than in reality. 

There are also some common problems with the 002 (hot) 
profiles (Figures 6-14, 6-16, 6-18, 6-20, 6-22, 6-24). 
Observations show an elevated inversion layer at 05/00Z and 
06/002, evidence of cooler marine air advancing inland. 

The forecasts in each experiment in each domain fail to 
recognize this feature. 

Model forecasts show winds strongly backing between 
layers two and four to a definitively offshore direction in 
layers five and six. These layers are also where model 
wind speed peaks for each experiment except COLD. 

Observations for 05/00Z and 06/00Z show peak speed in 
lower layers and backing much less with height. This is a 
typical signature expected for a strong sea breeze regime, 
very different from the signature presented by any of the 
forecasts. 

Lastly, the 3km profiles are very close to the 9km 
profiles, similar to the nighttime comparisons. Neither 
domain shows a greater propensity than the other to match 


the observed profiles. 
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VII. DISCUSSION AND RECOMMENDATIONS 


A. COMMON WEAKNESSES 

Tables 6-1 through 6-4 present an important 
consideration for this study: no model domain, in any 
experiment, forecasts with distinguishable skill. No 
forecast wind speed, wind direction, or temperature error 
is less than the standard deviation of the observations. 
To show skill, numerical model error should at least be 
less than the standard deviation of the observations 
(Kuypers, 2000). 

Another weakness seen throughout the study is the 
difficulty that each domain forecast has in approaching the 
cold part of the diurnal cycle. Forecasts tend to show 
higher winds at night, a result of the inaccurate lower 
layer thermal structure. Stronger winds would produce a 
deeper mixed layer and less pronounced nighttime inversion. 
Accurately depicting and tracing the evolution of the 
marine boundary layer is a critical element of pollutant 


forecasting (Ulrickson and Mass, 1990). 


B. HORIZONTAL TEMPERATURE STRUCTURE 

Figures 7-1, 7-2, and 7-3 show the nighttime CTRL 
near-surface horizontal temperature structure 12 hour 
forecasts over the entire 3km domain. Note the large areas 
of excessively warm forecasted temperatures shoreward of 
the basin's surrounding mountains, in some places more than 
10C greater than the observed temperatures. 

Figures 7-4 through 7-6 show the 22m AGL forecasted 
winds for the same times. Note on these figures the areas 


of downslope flow near the surface collocated with the warm 
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pockets seen in Figures 7-1 through 7-3. Subsidence 
associated with downslope flow acts to enhance the pre- 
existing coastal temperature gradient. This subsidence 
causes warming, due to adiabatic compression, which occurs 
in the warm air over land. 

The strong coastal surface temperature gradient 
produces a strong thermal wind. Lower layer geostrophic 
wind plus thermal wind equals upper layer geostrophic wind. 
This large thermal wind is what drives the stronger upper 
wind speeds in CTRL than what was actually observed (Figure 
6-15). The thermal wind is also what makes the model flow 
direction aloft coast parallel compared to onshore flow in 
the observations. The same can be said for forecasts in 
TEST1 and TEST2. COLD shows the least difference from 
observed temperatures in the lower layers. COLD, 
therefore, would have a less pronounced coastal temperature 
gradient, resulting in a weaker thermal wind. We 
hypothesize that this is why COLD's vertical speed 
structure more closely matches observations at 122. 

Figures 7-7a, 7-7b, and 7-7c show the 3km domain 12 
hour vertical profile forecasts for profiler site VNS. 
Tables 7-1 and 7-2 give the actual temperature and wind 
measurements taken from the VNS profiler at these same 
times, as well as some key lower-level model forecast wind 
values on Table 7-2. Note on Table 7-1 the slight increase 
in temperatures over the 48 hour period, again showing the 
general warming trend mentioned earlier. Table 7-2 shows 
the large disagreement between observed and forecasted 
winds at VNS. Downslope forecasted wind speeds near the 
surface are 253 times greater in the model than che 
observed downslope winds at VNS. 
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The model forecasted Skew-T plots do a fair job in 
capturing the shape of the temperature lapse when compared 
against the profiler temperature observations shown in 
Table 7-1. The forecasts show a profile indicative of 
subsidence at the surface, and the observations 
substantiate this with inversions between 473m and 683m AGL 
during the three 12Z times. 

Figures 7-8, through 7-10 show a N-S vertical cross 
section through the VNS profiler site, from north of the 
site down south to the coast (section depicted on Figures 
7-4 through 7-6). These figures again depict the 3km 
domain 12 hour forecasted temperatures and windspeeds. 

At 04/122 (Figure 7-8), the 31C isotherm sits aloft 
over VNS. The layer of temperatures greater than 31C over 
VNS increases each successive night (Figures 7-9 and 7-10). 
There is also a slight change in the 700mb wind forecasts 
over the three 12Z times. Higher heights further inland 
and lower heights offshore shift their orientation with the 
respect to the coast. The winds at 700mb back from 
slightly offshore at 04/12Z (Figure 7-11) to coast-parallel 
at 05/12Z (Figure 7-12) to slightly onshore at 06/12Z 
(Figure 7-13), reflecting a changing synoptic-scale 
pressure gradient orientation relative to the coastline. 

The warm low-level temperatures near the coast at 
04/122 (Figure 7-8) contribute to decreased surface 
pressure near the coast by effectively reducing the mass of 
the surface-to-700mb layer. This results in an increased 
pressure gradient along the coast. The increased gradient 
results in increased downslope winds toward the coast. 

This leads to speed convergence near the coast, and 


increased upward vertical motion due to continuity. The 
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increased downslope flow also produces more warming over 
areas near the coast due to adiabatic compression. This 
warmer air will contribute to decreased surface pressure 
near the coast, and the self-induced feedback loop will 
continue. 

The motion arrows in Figure 7-8 depict circulation 
vectors in the cross section. The 12 hour forecast valid 
at 04/122 shows near surface windspeeds in excess of 3 ms'^! 
near VNS, far greater than indicated by the observations in 
Table 7-2. Figures 7-9 and 7-10 show 12 hour forecasts 
valid at 05/122 and 06/122. At 05/122, the increased near- 
surface convergence produces greater upward vertical motion 
and is collocated with a low-level jet (LLJ). When the 
convergence decreases, as in Figure 7-10, the LLJ also 
decreases. The peak downslope winds, observed and 
forecasted at 05/12Z, correspond with maximum near-surface 
convergence and maximum strength of the LLJ. 

The maximum strength of the LLJ occurs when the 
thermal gradient is parallel to the synoptic geopotential 
gradient. Thus the synoptic fields (Figure 7-12) and 
mesoscale pressure fields are working constructively, 
resulting in greater wind speeds in the lower levels (Table 
7-2). 

A preliminary look at long-term forecasts (15-24hr) 


shows similar subsidence problems for 24 hour forecasts 


valid at 12Z. 


Се FINER RESOLUTION FORECASTS 
An important consideration in any model is the 
significance of forecast resolution and resolution's 


influence on a forecast. Finer resolution models need 
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finer scale observations in order to exploit the forecast 
strengths built into the model. 

For this study period, there is a definite gain in 
using the 3km or 9km domain vice the 27km domain. From 
Tables 6-1 through 6-12, no 27km forecast statistics are 
better than both the 3km and 9km forecasts. However, it is 
difficult to find any recognizable gain using a 3km domain 
forecast vice a 9km domain. The statistical results 
suggest little difference between the 9km and 3km 
forecasts. In each experiment, the 3km and 9km domain 


vertical structure plots are nearly identical. 


De MODEL PHYSICS 

Each experiment is a change, to some degree, of the 
model physics. Although there was no experiment that 
forecasted particularly well, COLD presented the "least 
bad" experiment. For COLD to outperform CTRL, the pre- 
forecast had to inject more error into the model 
initialization. This error presents a cascading effect. 
The forecast has error, is used in the next time step's 
initialization, and the model has great difficulty in 
correcting itself. Theory says COLD loses any mesoscale 
structure continuity from model run to model run, but this 
study shows it also minimizes perpetuating mesoscale error. 

The less complex PBL scheme in TEST1 was outperformed 
by the more sophisticated physics of CTRL. This study 
substantiates the theory that more complex physics schemes, 
when based on correct observations and valid assumptions, 
better resolve and forecast mesoscale structure. 

Adding observations did little to improve forecast 


results. Inserting extra surface observations into the 


Bo 


model initialization should give a more complete surface 
picture to the model. However, it does little to change 
the vertical structure in the initialization process. The 
low-level vertical structure is the key to accurate PBL 


analyses and forecasts. 


E. RECOMMENDATIONS 
This study presents interesting results, but poses 
many more questions worthy of study. Some recommendations 


for future studies: 


т; Try running MM5 using different soil models. 
will a different soil parameterization scheme result 
in better surface temperature forecasts valid at 122? Will 


it affect the coastal temperature gradient? 


2& Try initializing MM5 with profiler observations. 
Will adding upper air observations give more accurate 
initial vertical structure? Will better vertical structure 


at analysis time produce better forecasts? 


Sy Increase the number of vertical levels in the 
model boundary layer. 
Will greater boundary layer definition improve the 


performance of different PBL schemes and MM5? 


4. Apply these experiments at a different time of 
year. 
Do we have model problems that are linked to a 
specific synoptic-scale pattern? Are there similar model 
problems under different seasonal or synoptic conditions? 
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5. Employ a different initialization scheme (3DMQ, 
3DVAR, 4DVAR). 
Does analyzing fields three-dimensionally better 


capture and maintain realistic mesoscale structure? 


B 
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Table 3-1. Observation stations. 
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Table 4-1. Division of atmospheric layers. 
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Tables 6-1a, 6-1b, and 6-1c. CTRL case statistics. 
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Tables 6-2a, 6-2b, and 6-2c. COLD case statistics. 
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TESTI Case 
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Tables 6-3a, 6-3b, and 6-3c.  TEST1 case statistics. 
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TEST2 Case 
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Tables 6-4a, 6-4b, and 6-4c. TEST2 case statistics. 
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m AGL Temp 
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Table 7-1.  VNS profiler virtual temperatures for 04/12Z, 05/122, and 
06/122. 
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$ 34.22М 118.49W 


04/122 os/122 ^ [96/122 
m AGL| Wspeed (m/s) Wspeed (m/s) WindDir | 
| oss | ms | ons | ms | ons | ms | ons | ms | ons | ms | ons | ms | 
152 | 0.73 
| 210 | 0.95 | 3.00 | 18.77 | 47.80 | 2.47 | 6.70 |346.05|48.40| 1.73 | 2.40 | 55.96 |113.00| 
[326 | 0.76 | 4.70 | 14.23 | 88.50 | 5.15 | 9.10 |348.79|65.60| 2.77 |4.00 | 44.34 | 90.40 | 
84 1.48 34.24 | 99.00 2265 349933] 70.40] 2.97 4.50 41.52 | 88.50 | 
| 
| 
| 558 | 3.01 | 8.10 |113.96|109.50| 7.20 | 8.50 | 2.24 |78.80| 5.77 | 5.80| 5.97 | 82.30. 
790 | 2.38 | [252.90] | 3.24 6.90 | 11.54 | 81.80| 4.97 9.12 | 81.50 
848 3.49 251.73 3.66 2766 5.12 15358 81,70 


| 906 | 2.54 | _ [244.73| | 3.89 | [0.78 | — [4.90 | 6.60 | 19.39 | 81.90. 
[3.75 | — (|241.74| — | 4.04 | _ [350.78 [4.21 | 6.60 | 351.62 | 82.10 |. 
1022 | 4.26 245.93| | 4.28 | 337.57 4.83 |6.70 | 10.66 | 82.40 № 
| 1080 | 3.86 | — |248.80| || 3.71 | |мл) (| 5-46 | ТЕБЕ 
[1138 | s.os |  j2s2.as]_ | x.64 | — [323.34| [4.48 | — |ss0.38| — | 
1196 | 5.68 | |256.75| [4.20 | _ |зі2.64| | 5.48 | —— | 13.38- 
[1254 |999.00| ^ |9990) | 3.27 | [308.10]  |5.69 | — [359.17 
[1312 | 4.86 | — [242.08] ἰ3.88]᾽ |2926) [6.63 | [0.84 | | 
[2370 | 4.67 | [237.57] | 3.60 | —— 1277.18| —— | 5.90 | ]348.54᾽ 
1428 | 4.24 | _ [173.65] _ | 3.24 | — ]2z1.82| ^ | 6.24 || ОВ И 
1486 | 4.41] [173.94] [3.75 | [273.59] | 4.2013 | | 6.58 | — | 
[1544 | 4.53 | — [177.63| | 2.87 | |262.13| | 2.95 | — [|356.933] 2, 
[1602 | 3.91 | — [174.53] _ | 2.35 | [257.321 [2.31 | — | x38 | . .. 
[1660 | 4.14 | — |177.05| _ | 2.89 |__ [251.80 — [2.65 | [96 | 
1728 | 3.80 |  |178.76| 1.78 240.30 | 3.26 | 135 
[1776 | 4-01 | —— [|196.47| —— | 2.42 | [216.75] |36) EEEN И 
[1834 | 3-33 | —— [187.17| — (| 1.01 | [202.52| ο ἵ-- 
[1892 | 3.37 | — |189.00| — | 1.65 | —— [182.46| — | 2.78 | — |21.09 |] — — 
1950 | 2.94 189.45 1.69 172.06 ESL 31.59 
[2008 | 2.94 | — [179.13| — | 2.09 | [166.309] ος]. [asas] ᾿ 
[2066 | 2.85 | — |179.69! — | 1.97 | —— [162.79| 11.73] ἠ204.43᾽ 
[2124 | 3.22 | — lisa. sel [|2.55 | — Ira) | 3.88 | — [|220.08| — | 
3.64 188.35 146.62 2.44 
[2240 | 3.09 | _[191.42! _[ (| 2.45 [166260 [iz II 
2298 | 4.02 |192.]᾽ >50] hasa m= =u 


Table 7-2. VNS profiler winds for 04/12Z, 05/12Z, and 06/122 
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Figure 1-1. Numerical Model Forecast Error. 
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Figure 2-1.  Concentric domain boxes: 
3km (innermost). 
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81km (outermost), 





27km, 9km, and 
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Figure 3-1. Major topography of LA basin. 
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Figure 4-1. Mean surface temperature, 4 domains and observations. 
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Figure 5-1. Land Breeze / Sea Breeze schematic. 
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Figure 5-2. 04/002 NOGAPS 500mb heights(solid) and vorticity(dashed). 
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Figure 5-3.  04/12Z NOGAPS 500mb heights (solid) and vorticity (dashed). 
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05/122 NOGAPS 500mb heights(solid) and vorticity (dashed). 
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Figure 5-6. 06/002 NOGAPS 500mb heights (solid) and vorticity(dashed). 
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Figure 5-9. 07/122 NOGAPS 500mb heights (solid) and vorticity (dashed) 
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Figure 5-10. 08/002 NOGAPS 500mb heights(solid) and vorticity(dashed). 









Мон Апа 04 2997 90007000 NOG_ldeg (925 pres (C)* 
MON Au 1997 007000 Нос 14е4-<9 pres) Geopotenti Height 


Figure 5-11. 04/00Z NOGAPS 925mb heights (dashed) and temps (solid). 
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Figure 5-15. 06/002 МОСАРЅ 925mb heights (dashed) and temps (solid). 
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Figure 5-16. 06/122 NOGAPS 925mb heights(dashed) and temps (solid). 
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Figure 5-19. 08/002 NOGAPS 925mb heights(dashed) and temps (solid). 
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Figure 5-20. 04/002 NOGAPS surface mean sea level pressure. 
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Figure 5-21. 04/12Z NOGAPS surface mean sea level pressure. 
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Figure 5-22. 05/002 NOGAPS surface mean sea level pressure. 
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Figure 5-23. 05/122 NOGAPS surface mean sea level pressure. 


63 





1012 
1020 
Ly 


186 QS igi 
1452 NIE 


1020 


МЕр Апд 06 199%-0000У000 NOG ldeg (SFCj Sea Leyél pressura —.. 


Figure 5-24. 06/002 NOGAPS surface mean sea level pressure. 
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Figure 5-25. 06/12Z NOGAPS surface mean sea level pressure. 
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Figure 5-26. 07/002 NOGAPS surface mean sea level pressure. 
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Figure 5-27. 07/12Z NOGAPS surface mean sea level pressure. 
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Figure 5-28. 08/002 NOGAPS surface mean sea level pressure. 
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Figure 5-29. Mean temperature, observations. 
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Figure 5-30. Mean wind speed, observations. 
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Figure 5-31. Mean wind direction, observations. 
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Figure 6-1.  CTRL mean wind speed. 
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Figure 6-2.  CTRL mean wind direction. 
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Figure 6-3.  CTRL mean temperature. 


Wind Speed Mean, Sfc Layer COLD 


(m/s) 





VOS o9 ab BD 5% . de 49 HG ον «qo «; „а, % nh «bb о «4 < 
^» (o $ fy 5 % ^v Q e ο 9 η, ° % Q => τον 9 η, ° % ^v Q 
Q Q Q mS y “м η, су Q Q O Sy “м “у T © © Q N N sx η, © 
% > > > > > > ee < <? ee < < e e < «Ὁ ο ce < ce e < e 
Time 
——B— MM5 3km —— MHS 9km —s— YOCSS —e— ENS 27km —&— HNS 8ikm - 20» -Obs 


Figure 6-4. COLD mean wind speed. 
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Figure 6-5. COLD mean wind direction. 
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Figure 6-6. COLD mean temperature. 
70 


Wind Speed Mean, Sfc Layer TEST1 


AAA Edi 
A E A 

ANAL BA NET IM. "ND 
TAA ОАА ПО У 
EXA A MIE ТИС МА 
ТТНТ КЫ 


ES 


— 


PFF IE Ут е 


Тіте 





Figure 6-7. TEST1 mean wind speed. 
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Figure 6-8.  TEST1 mean wind direction. 
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Figure 6-9. TEST] mean temperature. 
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Figure 6-10.  TEST2 mean wind speed. 
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Figure 6-11.  TEST2 mean wind direction. 
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Figure 6-12.  TEST2 mean temperature. 
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Figure 6-13. 04/122 9km mean vertical profiles. 
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Figure 6-14. 05/002 9km mean vertical profiles. 
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Figure 6-15. 05/12Z 9km mean vertical profiles. 
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Figure 6-16. 06/002 9km mean vertical profiles. 
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Figure 6-17. 06/122 9km mean vertical profiles. 
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Figure 6-18. 07/002 9km mean vertical profiles. 
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Figure 6-19. 04/12Z 3km mean vertical profiles. 
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Figure 6-20. 05/002 3km mean vertical profiles. 
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Figure 6-21. 05/12Z 3km mean vertical profiles. 
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Figure 6-22. 06/00Z 3km mean vertical profiles. 
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Figure 6-23. 06/122 3km mean vertical profiles. 
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Figure 6-24. 
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Figure 7-1. 04/12Z MM5 CTRL 3km 12 hour temperature forecast for sigma 
level 9970 (approx. 22m AGL) with surface observations. 
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05/122 MM5 CTRL 3km 12 hour temperature forecast for sigma 


level 9970 (approx. 22m AGL) with surface observations. 


Figure 7-2. 
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Figure 7-3. 06/122 ММ5 CTRL 3km 12 hour temperature forecast for sigma 
level 9970 (approx. 22m AGL) with surface observations. 
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Figure 7-4. 04/122 MM5 CTRL 3km 12 hour wind (kts) forecast for sigma 
level 9970 (approx. 22m AGL). 
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Figure 7-5. 05/122 MM5 CTRL 3km 12 hour wind (kts) forecast for sigma 
level 9970 (approx. 22m AGL). 
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Figure 7-6. 06/127 MM5 CTRL 3km 12 hour wind (kts) forecast for sigma 
level 9970 (approx. 22m AGL). 
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Figures 7-7a, 7-7b, and 7-7c. MM5 CTRL 3km domain vertical profile 
12 hour forecasts at same lat/long of VNS profiler. 
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04/122 MM5 CTRL 3km N-S cross section forecast. 
Direction (arrow) and scalar speeds 


Temperature in Celisus (dashed). 
Isotachs (solid) are positive (southward) or negative 


(northward). 


Figure 7-8. 
in m/s. 
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06/122 ММ5 CTRL 3km N-S cross section forecast. 


Temperature in Celisus (dashed). 


Figure 7-10. 
in m/s. 


Direction (arrow) and scalar speeds 


Isotachs (solid) are positive (southward) or negative 


(northward). 


LLJ (heavy solid) isotach 10 m/s (out of page). 
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Figure 7-11. 04/122 MM5 CTRL 3km 12 hour 700mb height(m) and wind(kt) 
forecasts. 
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Figure 7-13. 
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